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Interacting quintessence models have been proposed to explain or, at least, alleviate the coinci- 
OO ■ dence problem of late cosmic acceleration. In this paper we are concerned with two aspects of these 

' kind of models: (i) the dynamical evolution of the model of Chimento et al. [L.P. Chimento, A.S. 

O \ Jakubi, D. Pavon, and W. Zimdahl, Phys. Rev. D 67, 0835f3 (2003).], i.e., whether its cosmolog- 

, ical evolution gives rise to a right sequence of radiation, dark matter and dark energy dominated 

eras, and (ii) whether the dark matter dark energy ratio asymptotically evolves towards a non-zero 
' constant. After showing that the model correctly reproduces these eras, we correlate three data sets 

that constrain the interaction at three redshift epochs: z < 10'', z = 10'', and z — 1. We discuss the 
model selection and argue that even if the model under consideration fulfills both requirements, it 
is heavily constrained by observation. The prospects that the coincidence problem can be explained 
by the coupling of dark matter to dark energy are not clearly favored by the data. 
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^ ' Recent measurements of luminosity distances using supernovae type la (SNIa) pj], of the cosmic microwave 

background (CMB) temperature anisotropies with the Wilkinson Microwave Anisotropy (WMAP) satellite 0], 
'nI" ' large scale structure "s], the integrated Sachs- Wolfe effect [3], and weak lensing fs], strongly suggest that 

^ , the Universe is currently undergoing a phase of accelerated expansion -see 6] for recent reviews. Within 

general relativity the obvious candidate to explain the present acceleration is the cosmological constant (or 
vacuum energy), and in fact the na'ive model built on it, ACDM, seems to pass reasonably well all cosmological 
tests. However, it suffers from two serious drawbacks from the theoretical side: the unnatural low value of 
the corresponding energy density, 123 magnitude orders larger than observed, and the so-called "coincidence 
problem", namely, "Why are the densities of matter and vacuum of the same order precisely today?", that 
' requires the vacuum energy density to be 96 orders of magnitude smaller than the matter density at the Planck 

, scale. (It is fair, however, to mention the existence of proposals in which a vacuum energy density of about the 

right order of magnitude stems from the Casimir effect at cosmic scales -see \7\ and references therein). This is 
why much attention has been devoted to models featuring an evolving and nearly un-clustered form of energy, 
■ usually dubbed "dark energy " , possessing a strong negative pressure high enough to drive late acceleration -see 

H ' [1| for an ample review of models. 

For simplicity, most cosmological models assume that matter and dark energy interact gravitationally only. In 
the absence of an underlying symmetry that would suppress a matter-dark energy coupling (or interaction) there 
is no a priori reason to dismiss it. Further, the coupling is not only likely but inevitable Q and its introduction 
is not more arbitrary than assuming it to vanish. On the other hand, it may help explain the coincidence 
problem. Ultimately, this question will have to be resolved observationally. Among other things, the interaction 
can push the beginning of the era of accelerated expansion to higher redshifts and it may erroneously suggest 
(if the interaction is ignored when interpreting the data) an equation of state for the dark energy of phantom 
type -see jTo| and references therein. Another question, is the form of the coupling. There is no clear consensus 
on this point and different versions, that arise from a variety of motivations, coexist in the literature. 

Cosmological models where dark matter (DM) and dark energy (DE) do not evolve separately but interact 
with each other were first introduced to justify the small current value of the cosmological constant dl!| and 
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nowadays there is a growing body of literature on the subject -see, e.g. [T4I and references therein. Recently, 
various proposals at the fundamental level, including field Lagrangians, have been advanced to account for the 
coupling [13| . Lagrangians of scalar field coupled to matter generically do not generate scaling solutions with a 
dark matter dominated period lasting long enough as required by cosmic structure formation 

In this paper we compare the interacting quintessence model (IQM) of Chimento et al. (see also its 

forerunner [3]) with observational data (supernovae, cosmic microwave background (CMB), and matter power 
spectrum) to set limits on the strength of the interaction DM/DE. The model was built to simultaneously 
account for the late phase of acceleration in the framework of Einstein relativity and significantly alleviates the 



coincidence problem. It evades the limits set in |14| and is compatible with a right succession of cosmic eras - 
radiation, dark matter, and dark energy dominated expansions. In a recent paper, Guo et al. also set constraints 
on interacting quintessence models However, the interactions studied by these authors differ from the one 
in [l5| . and while they use the cosmic background shift parameter and baryon acoustic oscillations alongside 
supernovae data to constrain the interaction, they do not consider the matter power spectrum, whereby our 
analysis may be viewed as complementary to theirs. 

The outline of this paper is as follows. Next section studies the critical points of the autonomous system of 
equations associated to the IQM [l^. Section III considers the restrictions set by Amendola et al. on the 
model to conclude that the latter evades these restrictions and, in particular, that an early epoch of baryon 
dominance is possible only if the strength of interaction is unnaturally large (beyond the limits set by the CMB 
data). Section IV focus on the observational bounds coming from the CMB, matter power spectrum and recent 
supernovae type la data. Finally, section V summarizes our results. 



II. DYNAMICS OF THE INTERACTING QUINTESSENCE MODEL 



If the quintessence DE decays into DM, both energy densities evolve differently than in non-interacting 
quintessence cosmologies and, therefore, the interaction can be tested by its effects on the dynamical evo- 
lution of the Universe. Due to the interaction, the fraction of cold dark matter (CDM) at any time in the 
past is smaller than in non-interacting models with the same cosmological parameters today. Since the dark 
matter energy density grows more slowly the beginning of the period of accelerated expansion and the temporal 
evolution of the gravitational potential differ from non-interacting models. Observables such as the angular and 
luminosity distances depend on the time evolution of the energy density. But the effect does not only occurs 
at zeroth order; the evolution of first order matter density perturbations is also affected and so is the pattern 
of anisotropics of the CMB. This section describes the dynamical evolution of dark matter and dark energy 
densities at zeroth order in the interacting quintessence model of Ref. [15] to point out the main differences 
with respect to models with no interaction. 

The phenomenological model of Chimento et al. (see also [iB]) assumes that the dark energy decays into 
cold dark matter thereby both energy densities do not longer evolve separately. In [l^ the following ansatz 
was adopted 

Pc + iHpc = 'iH c^{p^, + Pc) , px + iH(l + Wj;)px = -3H c^(p^, + pc), (1) 

where p denotes the energy densities (subscript c for cold dark matter and subscript x for dark energy), Wx is 
the equation of state parameter of dark energy, H = a^^da/dt the Hubble function and a the scale factor of 
the flat Friedmann-Robertson- Walker metric. The model assumes that the dark energy decays just into cold 
dark matter and not into any other component such as neutrinos, baryons, or photons. Decay into neutrinos 
have been recently studied in the literature p^ . The coupling with baryons is constrained by measurements 
of local gravity [13, [HI . The above ansatz is essentially phenomenological; it is the simplest interaction model 
that leads to a constant ratio at early and late times which certainly alleviates the coincidence problem. While 
at present a rigorous derivation from first principles does not exist, the ansatz may be roughly justified as 
follows: the right hand side of both equations in ^ must be functions, say Q and —Q, of the energy densities 
multiplied by a quantity with units of inverse of time. For the latter the obvious choice is the Hubble factor 
H, so we have that Q = Q{Hpx ,Hpc). By power law expanding Q and retaining just the first term we get 
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Q ~ Aa; Hpx + XcHpc- To facilitate comparison of the model with observation it is suitable to eliminate one 
the two A parameters. Thus we set A^; = Ac = 3c^ and arrive to Eq. ([1]). The simpler choice \c — would not 
yield a constant dark matter to dark energy ratio at late times. The dimensionless term 3c^ measures to what 
extent the decay rate differs from the expansion rate of the Universe. Thus, the model is characterized by this 
single parameter, c^, which also gauges the intensity of the interaction. The lower c^, the closer the evolution 
of the Universe to a non-interacting model is. 



Equations in ([T]), the continuity equations for photons and baryons and Friedmann equation = {k^ /'S){prad + 
Pb + Pc + Px), constitute a closed system. In Fig. [T] we plot the evolution of the different matter density 
components: cold dark matter (thick solid line), dark energy (thin solid), photons (thin dot-dashed) and baryons 
(thick dot-dashed) for two values of the interaction parameter: left and right panels correspond to — 10^^ and 
= 0.1, respectively. In the figure we can see that Eqs. ([T]) have two asymptotic solutions with r = constant: 
at early times and in the immediate future. It is easy to prove that, in these two cases, -I- (w^-c"^ -|- 2)r -1-1 = 
and the constant ratio is r_ — r{z — )■ oo) ^ c^^ and r+ = r(z 0) ^ c^. These results represent a clear 
conceptual advantage with respect to the ACDM model. If the present acceleration of the Universe is generated 
by a cosmological constant, then the initial values of dark matter and cosmological constant have to be tuned by 
96 orders of magnitude at Planck time, i.e., r{tpianck) = ^m/^K ~ 0(10^^). In the model described above this 
ratio is fixed by the interaction rate, c^. The initial condition problem is significantly alleviated. For example, 
if '-^ 10""* the initial dark matter to dark energy ratio will be r{tpianck) ~ 0(10**). 



To study the evolution of the corresponding autonomous system formed by the continuity equations of all 
energy components and Friedmann's equation we introduce the following set of variables. 



H Px K Pc K Pb K prad 



(2) 



ith K = vSttG. Then, the dynamical equations can be recast as 



Zwx [x^ -l)+u^ -Zc^[l + —. 



X, 



SwxX-' + + Sc" I 1 H — - 



- [iWxX^ +U^]z , 



(3) 



= [iWxX^ + - l] U , 



x"^ + y'^ + z"^ + V? . 



where a prime denotes derivative with respect to In a. As said above, the motivation for interacting quintessence 
models is to solve or, at least, ameliorate the coincidence problem. A solution will be achieved if the system 
(121) presents scaling solutions [23]. As already discussed, scaling solutions are characterized by a constant dark 
matter to dark energy ratio. Even more important are those scaling solutions that are also an attractor. In 
this way, the coincidence problem gets substantially alleviated because, regardless of the initial conditions, the 
system evolves toward a final state where the ratio of dark matter to dark energy stays constant. 

In Table U we give the location of the four critical points of the autonomous system of Eqs. ([21), their 
stability character and whether they give rise to accelerated or decelerated expansions. The first critical point, 
(a;?^, yL, 0, 0), corresponds to the dark energy dominated epoch with 

xl_ = [wx - Ac?Wx^ , 2/^ = 1- (wx - 4c2u;:r) , (4) 

and it requires < |it;a;|/4. As shown in p^ . this coincides with the condition the interaction must fulfill to 
give a physically acceptable evolution for the Universe. 
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Existence condition 


Stability character 


Acceleration 


{xl,yl,0,0) 


< < \w^\/4 


attractor 


a > il -U^a: < — 1/3 


(4,y^,0,0) 


< < \w^\/4: 


saddle point 


a < Q\/c^,W:x 


(0,0,1,0) 




unstable 


a < 


(0,0,0,1) 




unstable 


a < 



TABLE I: Location of the critical points of the autonomous system of Eqs. ([3]), their stability and dynamical behavior 
of the Universe at those points. 

The second critical point, 



= - — [wx + \/wx+ '^c^Wx] , y\ = l- - — [wx + \/wl+ '^C^Wx] , (5) 

IWx ^ ' iWx ^ ' 

corresponds to the CDM dominated era and it is a saddle point. 

The third critical point, (0,0,1,0), is unstable and physically unrealistic. It corresponds to a universe with 
just baryons (it contains neither radiation, nor matter, nor dark energy). Finally, the fourth point, (0,0,0,1), 
corresponds to the radiation dominated era and it is also unstable. To summarize, as soon as the other energy 
densities grow more important than radiation, the system moves away from this point, reaches a CDM dominated 
era, which corresponds to a saddle point, and eventually moves to a final attractor scaling solution in the dark 
energy dominated era. In a natural way, the dynamical evolution of the Universe is to asymptotically approach 
this attractor, a never-ending phase of accelerated expansion, in which the energy ratio r = pd Px remains 
constant. 



III. IQM AND THE COINCIDENCE PROBLEM 



In jl4| Amendola et al. discussed the conditions that interacting quintessence models must verify in order 
to solve the coincidence problem and have a correct sequence of cosmological eras -i.e., radiation, matter, and 
accelerated scaling solution. Many of the Lagrangians proposed in the literature do not correspond to acceptable 
scaling solutions. Those models either enter a phase of accelerated expansion right after the radiation epoch or 
present baryon dominated periods that would affect the growth of structure in the Universe. This brings about 
a great stress on the ability of a large class of models to explain or even ameliorate the coincidence problem. 

The phenomenological model of Chimento et al. [l^ is not affected by the results of Ref. [l3| since the coupling 
between dark matter and dark energy is not of the type studied by Amendola et al., namely, Q cx pm d(j) / d{hi a) , 
where pm denotes the dark matter energy density and (j) the scalar field. Equivalently, the Lagrangian associated 
to the effective potential given in Eq. (5) of Ref. [l^ is not of the type described in [l3|. As can be seen in Fig. 
[T]the model contains two scaling regime solutions, one in the radiation- and matter-dominated eras, another in 
the period of accelerated expansion -at late times. It reproduces the right succession of radiation, matter and 
accelerated epochs, each lasting long enough to give rise to the observed cosmic structure. It only remains to 
show under which circumstances there is not a baryon dominated period (as in left panel of Fig. [1]) so that the 
formation of galaxies and large scale structure go unaffected. 



A. Baryon dominated era 



Some interacting quintessence models present the severe drawback that the Universe undergoes a baryon- 
dominated period prior to the DM era. This is rather problematic since CMB temperature anisotro pics and 
the evolution of matter density perturbations would result largely affected. Unlike the models of Ref. [2J|, the 
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Chimento et al. model fisll does not have a stable baryon dominated era. To see this let us estimate when and 
for how long could the model we are considering undergoes a baryon dominated epoch. The maximum baryon 
energy density can be obtained by setting z' = in Eqs. ([3]). Noting that v? = —'iwx , we have from the last 
equation in ^ that 

^max = 1 - - (1 - iWx)x^ . (6) 

For r — constant, + {wx + 2)r + 1 = 0. To affect the growth of structure and the CMB temperature 
anisotropics, baryon domination should occur, if at all, well before the period of accelerated expansion, i.e., 
when 1 < r ~ r+ with r+ = (l/2)[-(u;^ c"^ + 2) + yj{wx c^^ + 2)2 - 4] ~ \wx\c~'^ ■ Since y"^ = r Eq. (O 
becomes 

^la. ^ 1 - (1 + 3c2)y2 . (7) 

Therefore, for an early baryon dominated era to have existed we must have z^^^ ~ 1 — (1 + Sc^)^^ > + . 
Thus, bearing in mind that <C the condition for the said epoch to exist is 

A baryon dominated epoch may have occurred for ~ 0(10"^) or larger. However, this does not correspond 
to any stable critical point of the autonomous system and -as we shall see below- such high values of are 
ruled out by the CMB data. Therefore, for values of about 10"^ or lower the model of Ref. [l^ is free of a 
long period of baryon dominance. 

To conclude, the succession of expansion eras in the model we are considering can be summarized as follows: 
after the initial radiation dominated era, the Universe might have gone through a very short baryon dominated 
period (if at all) followed by a DM dominated epoch (corresponding to an unstable critical point), to finally 
asymptotically approach a regime with a constant dark matter to dark energy ratio, at late times, in the 
accelerated expansion era. 



IV. COMPARISON WITH OBSERVATIONAL DATA 



The IQM of Ref. [TB| evades the restrictions imposed by [T3| and provides a viable model to explain the 
coincidence problem by means of cosmological scaling solutions. The actual existence of an interaction must 
be elucidated by contrasting the model with observations. In Refs. [l^ and [13] we studied the effect of the 
interaction on the temperature anisotropics of the CMB, on matter power spectrum and on the luminosity 
distance. This opens up the possibility of using observational data to constrain the strength of the interaction 
-i.e., the parameter-, at different redshifts. At present, data are available on the power spectrum of matter 
density perturbations extracted from galaxy catalogs (25l . [2^ [27j , on temperature anisotropics ranging from the 
largest scales down to ^ 20 arcmin ,2,] and on luminosity distances up to redshift z ~ 1 measured using SNIa 
as standard candles Each data set probes the DE/DM coupling at different periods: CMB anisotropics at 
z ~ 10"^, the slope of the matter power spectrum at z < 10'* and the luminosity distance test at z ~ 1 and 
below. We carried out independent statistical analysis for each data set to see what constraints did the data 
impose on the interaction at different epochs. 

The effect of the interaction on the various observables can be determined numerically only. Physically, the 
interaction renders the potential wells shallower in the matter dominated regime and matter density pertur- 
bations grow slower [23| ; the gravitational pull on baryons prior to recombination weakens and the interaction 
alters the relative height of the acoustic peaks of the CMB spectrum and, finally, the interaction also changes 
the time evolution of the different energy densities, i.e., it affects the luminosity distance of standard candles. 
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40<Ho< 100 km/s Mpc"^ 


< ^Ich^ < 1 


< nth^ < 1 


-l<w^< -0.5 


0.5 < 71^ < 1.5 


0.5 <As< 1.5 


0.1 < e"^^ < 1.5 


-9 < logc^ < -1 











TABLE II: Priors on different parameters for WMAP, SDSS data analysis. MCMCs are constrained to take values within 
those intervals. The amplitude of the matter power spectrum yl^ = 1 is the amplitude fixed by the COBE normalization. 

Figure [2] shows these effects. In Figure [2ja)-l2ljc) we plot the matter, radiation power spectra, and the lumi- 
nosity distance for three cosmological models with the parameters of WMAP Syr data concordance model and 
interaction parameter = 0, 10~^,0.1 corresponding to solid, dotted and long dashed lines, respectively. We 
also plot the data (see below) used in the statistical analysis. In Fig. EJa), the effect of the interaction is to 
shift the matter radiation equality and strongly damp the slope of the matter power spectrum. As a result, a 
large interaction erases the matter power spectrum at small scales. In the figure, power spectra are normalized 
to Cosmic Microwave Background/Differential Microwave Radiometers (COBE/DMR). In Fig. [UJb) the inter- 
action changes the relative height of the acoustic peaks. At low multipoles the amplitudes also differ, since the 
slower evolution of the potential at low redshift gives rise to an Integrated Sachs- Wolfe component of different 
amplitude. This effect is discussed in [s^. In Fig. ^c) we see that the effect of the interaction is not very 
significant, and only when is as high as ~ 0.1 there is a tiny difference with non-interacting models that can 
be perceived at z '--^ 2. 

To constrain the model with observations, for every given set of cosmological parameters we generate the power 
spectrum of matter density perturbations and temperature anisotropies of the CMB using a modified version 
of the CMBFAST code [2g|. Our cosmological models are described by the following parameters: dark matter, 
baryon and dark energy densities (Sic, = ^ — ^c~ ^b)i dark energy equation of the state parameter {wx), 

Hubble constant (i?o), amplitude and slope of the matter power spectrum at large scales {As, ns), interaction 
parameter (c^), and the epoch of reionization measured by the optical depth of the Universe to CMB photons 
(r). We do not consider gravitational waves, massive neutrinos, or universes with curved spatial sections. 
Effectively, our parameter space is 8-dimensional. A brute force analysis of this parameter space would prove 
computationally very expensive. Instead, we chose to explore the parameter space using a Monte Carlo Markov 
chain (hereafter, MCMC) method, where the parameter domain is more heavily explored close to the best fit 
values. This method speeds up the calculation and at the same time explores fully the likelihood function 
and allows to determine accurately the model parameters and its error bars. We implemented the method as 
described in Ref . [1^ which was first proposed by Gelman and Rubin . For each data set we run 4 different 
MCMCs initiated at different regions of the parameter space. We followed the same practice than the WMAP 
team: the chain automatically stopped when a given degree of convergence {R < 1.2 in the notation of [29j ) 
was attained and when, at least, 40,000 models were computed. 



A. Constraints on IQMs from WMAP Syr data 



We run 4 MCMCs comparing models with WMAP 3yr data, using the same likelihood code than the WMAP 
team which is available at the Lambda Archive [SS^I . Table |TT] shows the intervals where parameters are let to 
vary. These priors are conservative in the sense that they include all known measurements and we do not expect 
these limits to have any effect on our final results. Still, we first made a rough exploration of the parameter 
space to ensure that the maximum of the likelihood function remains far away from the surface of the parameter 
space volume, so the likelihood hypersurface is not affected by our priors. We stop the program when the chains 
have satisfied convergence and well-mixing criteria. 

Figure[3]presents the joint confidence intervals at 68%, 95% and 99.9% for pairs of parameters after marginalizing 
over the rest. For convenience, the axis uses a logarithmic scale and it has been cut at < 10^^. In the 
figure, contours run parallel to the vertical axis, meaning that the data has only statistical power to set an 
upper limit on the interaction parameter. At the Icr confidence level, < 2.3 x 10""^. Figure [3] shows the 
same joint confidence contours for other pairs of parameters. In Table IIIII we presents the best fit values of 
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c 


< 2.6 X iU 


< 4.0 X iU 




U.iUO_Q 


0.24l»;f3 




2.22 ±0.07 




Wx 


< -0.90 


< -0.6 


Ho 


7l.3t\i 




ris 


94+0 04 

'-'•^^-0.05 




T 


0.09 ± 0.04 





TABLE III: Mean values of cosmological parameters and their la errors obtained by fitting WMAP Syr and SDSS LRGs 
data. For Wx only upper limits at the la confidence level are given. Values not quoted are weakly or not constrained by 
the data. 



the cosmological parameters obtained from WMAP Syr data alongside their la error bars. Those intervals are 
in rather good agreement with WMAP Syr results: the values on ilch^ are slightly more uncertain and on Hq 
more precise. The small differences are to be expected due to our different sampling of the parameter space 
[sH . The main difference lies in Wx- For this parameter, we can only set upper limits. The reason is that we 
restricted our analysis to quintessence models, Wx < —1, and did not consider phantom models. 

Let us remark that the shape of the likelihood function is not altered by the introduction of a new degree 
of freedom. It is well known that the CMB temperature anisotropics are highly degenerated and different 
combination of parameters give rise to very similar radiation power spectra. One would expect that one extra 
parameter would introduce extra degeneracies and modify the overall shape of the likelihood function. This is 
not the case and it simply reflects the fact that the best fit parameters are rather insensitive to the interaction 
(see Fig. An obvious explanation is that WMAP data does not allow for large departures from non- 

interacting models. The interaction changes the depth of the potential wells at the last scattering surface but 
this effect cannot be counter-balanced by a suitable variation of other cosmological parameters. The maximum 
of the likelihood function is unaffected and the contours are not shifted. 



B. Constraints on IQMs from SDSS and 2dFRGS data 



The matter power spectrum has been recently derived from two different galaxy catalogs: 2-degree Field Red 
Galaxy Survey (2dFRGS) '25'|, that comprises a sample of more than 2 x 10^ galaxies with measured redshifts, 
and the Luminous Red Galaxy sample (LRGs) of the Sloan Digital Sky Survey (SDSS) [26j]. A larger study, using 
the SDSS Data Release 5, and comprising more than half a million galaxies showed similar results [l^. In [2^ . 
the power spectrum from the SDSS was computed producing uncorrelated minimum-variance measurements in 
20 /c-bands of the clustering power and its anisotropy due to redshift-distortions, with narrow and well-behaved 
windows functions in the range 0.01 h/Mpc < k < 0.2 ft./Mpc. We used these data to fit the theoretical power 
spectrum, P(fc), calculated with a modified CMBFAST software and thus constrain the parameters of the 
cosmological model. Table HIl shows the priors imposed on the parameters in our MCMCs runs. In addition, we 
allowed the bias parameter to vary in the range 0.5 < b < 2.5. We run 4 different MCMCs till the criteria of 
convergence and well mixing were met. 

Figure [5] shows the joint confidence contours of with other representative parameters. Like in Fig. [3l contours 
run almost parallel to the vertical axis. In Table IIIII we give the cosmological parameters measured from this 
data set. We derived an upper limit of < 4.8 x 10"'^ at la confidence level. This constraint is weaker than the 
one provided by WMAP 3rd year data, but applies to a different range in redshift. This bound sets an upper 
limit on the value of averaged over the whole length of the matter dominated epoch. Varying Hq or the total 
matter content changes the matter radiation equality and shifts the maximum of P{k) to larger scales, but it is 
around the maximum where the data have less statistical power. Matter and baryons change the amplitude and 
location of baryon acoustic oscillations but useful constraints can only be obtained by separating the physics 
of the oscillations from that governing the overall shape of the power spectrum [s^l loosing the information on 
the interaction. As a result, we obtained no useful constraints on cosmological parameters except on c^. We 
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also run a set of MCMCs for the 2dFGRS data, but while the values of the main cosmological parameters were 
in very good agreement with those of Table IIIIl the constraints were even poorer than those obtained from the 
SDSS data, so we do not present them here. 



C. Constraints from SNIa data 



At low redshifts, only the dark matter and dark energy densities contribute substantially to the dynamical 
evolution of the Universe. Since the interaction alters their relative ratio, it affects the actual distance to 
objects. Distance measurements, like the luminosity distance, d^, obtained from standard candles like the 
SNIa, are sensitive to the dark energy parameters, fix, Wx, and the strength of the interaction, c^. The 
likehhood function is given by [H, 

, 1 2n 2fTj o 2\ M-^{^J-the{zi■,Ho,Wx,^x,c'^)-^J■obs.iY /^x 

£cxexp(--x), X {Ho,Wx,^x,c ) ^ — ^ . (9) 

In this expression, fi = Slnd^ + 25 is the distance modulus, cr^ is the uncertainty on the supernova redshift (in 
units of the distance modulus) due to its intrinsic motion and the peculiar velocities of the parent galaxy and 
^?o6s i error in the determination of the distance moduli due to photometric uncertainties. Subindexes 

the and obs indicate theoretical and observed uncertainties, respectively. Since the intrinsic luminosity of SNIa 
is not known, the distance moduli are insensitive to Hq, only 3 parameters (fix, Wx, and c^) are needed to 
characterize a model. We carried out a brute force analysis sampling homogeneously the parameter space: 
(—9 < logc^ < —1, < < 1, —1 < Wx < —0.5); each interval was subdivided in 10 steps. Figure [6] shows the 
confidence intervals for vs fix after marginalizing over Wx- We remark that, in this case, the data is almost 
insensitive to c^, as we could have expected from Fig[2fc). 



D. Bayesian model selection 



To solve the coincidence problem a decay of DE into DM was postulated and, consequently, a new parameter, 
c^, was introduced to gauge the strength of the coupling. In this subsection we explore whether the existence 
of this coupling (i.e., a non- vanishing value) it is also justified by a statistical description of the data. The 
introduction of further parameters will lead to a better fit to the data. This must be balanced against the loss of 
predictive power of the theory. The Bayesian Information Criteria (defined as BIC ~ + k log N [sll , where 
= —2 log/!, k the number of free parameters, and N the number of data points), penalizes the inclusion 
of additional parameters to describe data of small size. For the WMAP Syr data x^(c'^ = 0) = 11252.06 and 
X^(best model) = 11251.85, so the interacting quintessence model for = 2.3 x 10"'^ yields a slightly better fit 
to data than the corresponding non-interacting model. However, the application of the Bayesian criterion gives 
B/C(best model) = 11269 and BIC{c^ = 0) = 11266, so the introduction of this parameter is not statistically 
encouraged. This results disagree with our analysis of WMAP lyr data ^is']. It was the unexpected consequence 
of our choice of priors: we fixed all models to agree with the COBE/DMR normalization. Since WMAP lyr 
concordance model preferred a lower normalization, it was heavily penalized compared with interacting models. 
This discrepancy shows the importance of priors in Bayesian analysis. As argued in [s^, model selection is not 
a substitute for parameter fitting; a lower BIC does not imply the data rejects interacting models, just makes 
them unlikely. However, they still have a theoretical advantage since they can help explain the coincidence 
problem, while non-interacting models cannot. 
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V. CONCLUSIONS 



Let us briefly point out that while the IQM model does not solve the coincidence problem in full (as it cannot 
predict the present ratio of dark matter to dark energy [s^]) it substantially alleviates it. Indeed, if we take 
initial conditions at the Planck scale, in a model without interaction the ratio r at that epoch is about 10^^. 
If ~ 10~^, a value compatible with WMAP upper limit then, asymptotically, r{t 3> to) ~ 10^, where to 
is the present time. Initially r{t <C to) ~ 10"'^, so, while in a ACDM model the initial conditions have to be 
fine-tuned by 96 magnitude orders to get the present accelerated expansion, in the IQM the difference in orders 
of magnitude is just 6. A further advantage is that the model is not limited by the results of Ref. since it 
successfully predicts that the Universe undergoes three successive cosmic eras of expansion, namely, radiation, 
matter and dark energy dominance. If the interaction parameter were ~ 0(10"^) or larger, it could undergo 
an era of baryon dominance, but this is not a stable solution and, in any case, values larger than 2.3 x 10"'^ 
do not appear compatible with CMB data. 

The existence of a late accelerated expansion of the Universe requires Wx < —1/3, as in noninteracting models. 
Whether or not the present value of the ratio r coincides with the final attractor value depends on the parameters 
and Wx- These are unknowns that can be constrained by observation. Our analysis shows that the CMB, 
matter power spectrum and luminosity distance data only set upper limits on the DM/DE coupling and, 
statistically, the introduction of this extra parameter is not favored by the data. At present, CMB data sets the 
tightest upper limit on the strength of the interaction (c^ < 2.3 x 10~^). As Fig. [7] reveals, this upper bound 
can be translated into an excluded region (shaded area) for the dark matter to dark energy ratio, r. Solid lines, 
from top to bottom, correspond to = 10^^, 10^"', 10^"^ and 2.3 x lO^'^ that bounds the shaded area. As the 
figure shows, the smaller c^, the longer it takes to reach the attractor solution. With WMAP 3yr constraint, the 
attractor solution will be reached only in 10^^ yr from now. Our upper limits on put pressure on this type 
of interacting models that, combined with the more generic results of [l3|, show the difficulties that interacting 
models encounter into explaining the coincidence problem. If the new WMAP 5yr data release sets an even 
stronger constraint on c^, the IQM will loose much of its appeal. 
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(1+z) (1+z) 

FIG. 1: Evolution of the different energy densities with redshift. Left panel: (? = 10~^, right panel: (? = 0.1. Thick 
and thin solid, and thick, and thin dot-dashed lines correspond to dark matter, dark energy, baryons and photons, 
respectively. 




FIG. 2: Matter, radiation power spectrum, and luminosity distance for a cosmological model with Q,x = 0.74 and 
Wx = —0.9 and three different interaction parameters <? = 0, 10~^,0.1, corresponding to solid, dotted, and long dashed 
lines, respectively. The cosmological parameters are those of the fiducial WMAP Syr data. 




FIG. 4: Joint confidence intervals at the same confidence levels as in Fig. [3] for pairs of parameters after marginalizing 
over the rest. 
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FIG. 5: Contours as in Fig. |3] obtained using the SDSS matter power spectrum data. 
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FIG. 6: Contours as in Fig. [3] obtained by fitting the luminosity distance IQM to the "Gold" sample of SNIa data of 
Riess et al. and SNLS data of Astier et al. 
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